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Abstract 

In the accompanied paper [14j, a delayed nonlinear model for pricing corpo- 
rate liabilities was developed. Using self-financed strategy and duplication 
we were able to derive two Random Partial DiflFerential Equations (RPDEs) 
r^ describing the evolution of debt and equity values of the corporate in the last 



delay period interval. In this paper, we provide numerical techniques to solve 
our delayed nonlinear model along with the corresponding RPDEs modeling 

O the debt and equity values of the corporate. 

^^ Using financial data from some firms, we compare numerical solutions 

from both our nonlinear model and classical Merton model |7| to the real 

• fh corporate data. From this comparison, it comes up that in corporate finance 

rS the past dependence of the firm value process may be an important feature 

C^ and therefore should not be ignored. 
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1. Introduction 

Due to the remarkable growth of the credit derivatives market, the inter- 
est in corporate claim value models and risk structure has recently increased. 
Financial distress tends to be an important factor in many corporate deci- 
sions. The two main sources of financial distress are corporate illiquidity and 
insolvency. In his paper [2j, Gryglewicz explains how changes in solvency 
affect liquidity and also how liquidity concerns affect solvency via capital 
structure choice. Corporate solvency is the ability to cover debt obliga- 
tions in the long run. Uncertainty about average future profitability, with 
financial leverage, generates solvency concerns. Corporate insolvency may 
lead to corporate reorganization or to bankruptcy of the firm in the worst 
case. Corporate bankruptcy is central to the theory of the firm. A firm is 
generally considered bankrupt when it cannot meet a current payment on 
a debt obligation. In this event, the equity holders lose all claims on the 
firm, and the remaining loss which is the difference between the face value 
of the fixed claims and the market value of the firm, is supported by the 
debt holders. In the literature of corporate finance, Merton |7j appears to be 
the main pioneers in the derivation of formulas for corporate claims. This 
model is a dual of Black and Scholes model [i2\ for stock price. Merton [7J 
further analyzed the risk structure of interest rates. More specifically, he 
found the relation between corporate bond spreads and government bond, 
and attempted to determine a valid measure of risk. He also developed the 
deterministc partial differential equation modelling the debt and equity of the 
firm. The assumption of constant volatility in the original Black-Scholes and 
Merton models from which most claims derivations are inspired, is incom- 
patible with derivatives prices observed in the market (see jSJ [U [6l [^ and the 
references therein). For stock price, two alternative theories are mostly used 
to overcome the constant volatility drawback. The first approach sometime 
called level-dependent volatility describes the stock price as a diffusion 
with level dependent volatility [lOj. The second approach sometime called 
stochastic volatility defines the volatility as an autonomous diffusion 
driven by a second Brownian motion. [^ 



^In the sense that the first Brownian motion drives the asset price 



In [T7j, a new class of nonconstant volatility model which can be extended 
to include the first of the above approaches, that we called delayed model 
is introduced and further study in P il6j for options prices. This model 
shows that the past dependence of the stock price process is an important 
feature and therefore should not be ignored. The main goal of this model 
is to make volatility self-reinforcing. Since the volatility is defined in terms 
of past behavior of the asset price, the self^reinforcing is high, precisely 
when there have been large movements in the recent past (see [17]). This is 
designed to refiect real-world perceptions of market volatility, particularly if 
practitioners are to compare historic volatility with implied. 

Following the duality between the stock price [12j and corporate finance 
[7] , we have recently introduced in [131 ISl the nonlinear delayed model in 
debt and guarantee. Using self-financed strategy and replication we estab- 
lished that debt value and equity value follow two similar Random Partial 
differential Equations (RPDEs) within the last delay period interval q The 
analytical solution of our nonlinear model and RPDEs are unknown in gen- 
eral case and therefore numerical techniques are needed. 

In recent years, the computational complexity of mathematical models 
employed in financial mathematics has witnessed a tremendous growth (see 
[221 [211 [T5] and references therein) . The aims of this paper is to solve nu- 
merically our delayed nonlinear model for firm market value along with the 
corresponding RPDEs, using real data from firms. Comparison will be done 
with classical Merton model. To the best of our knowledge such comparison 
has not yet been done in the financial literature. Two major comparisons 
will be preformed: the market value of each corporate and its equity value 
(or its debt value). We will first approximate the volatility of each corpo- 
rate, afterward solve numerically our nonlinear model for the market value 
of the corporate along with the corresponding Merton model using the 6— 
semi implicit Euler Maruyama scheme to obtain sample numerical solutions. 
Monte carlo method will be thereafter used to approximate the mean nu- 
merical solution of each model. The meam numerical value from each model 
(our nonlinear model and Merton model) will be therefore compared with the 
real market value {V) of the corporate. For debt value (or equity value) solu- 
tions of RPDEs established in the accompanied paper [14j, efficient numerical 
scheme based on finite volume-finite difference methods (discretization re- 



^The final time interval with length equal to the time delay 



spect to the firm value V) and exponential integrator (discretization respect 
to the time t) will be used. Recently, exponential integrators have been used 
efficiency in many applications in porous media flow [U EOj [IHl [25l [26] , but 
are not yet well spread in flnance. The same numerical technique is also used 
to solve deterministic Partial Differential Equations (PDEs) modeling debt 
value or equity value in Merton model. Comparisons are done with the real 
data from flrms for each model (our delay model and Merton model). 

From our comparison, it comes up that in corporate flnance the past 
dependence of the flrm value process is an important feature and therefore 
should not be ignored. The main goal of this paper is to call for further atten- 
tion into the possibility of modeling market value of the flrm with nonlinear 
delayed stochastic differential equations. 

The paper is organized as follows. In Section [!} we recall our delayed non- 
linear model for corporate claims as presented in [14j along with the Merton 
model [7j. In Section [3| numerical techniques for our delayed nonlinear model 
are provided. We flrst present the 9— semi implicit Euler Maruyama for the 
flrm market value V and provide numerical experimentations for both our 
nonlinear model and Merton model using real data for some flrms. We end 
this section by providing numerical technique to solve efflciently our (RPDEs) 
modeling the debt and equity of the flrm along with numerical experimen- 
tations for the two models (our delayed nonlinear model and Merton model) 
with real data for some flrms. The conclusion is provided in Section [4j 

2. Stochastic delay model for corporate claims 

Here we present the stochastic delay model formulated in the accompanied 
paper jT4] along with Random Partial Differential Equation (RPDE) that 
should satisfy any claim. We assume that: 

Ai The value of the company is unaffected by how it is flnanced (the capital 
structure irrelevance principle). 

A2 The market value of flrm at time t G [0,T], V{t)^ follows a nonlinear 
Stochastic Delay Differential Equation (SDDE) 

dV{t) = {aV{t)V{t -L)- C)dt + g{V{t - L))V{t)dW{t) 

(1) 
y(t) = (^(t), tG[-L,0] 



on a probability space (f^, J^, P)) with a filtration (^o<t<r) satisfying 
the usual conditions. 

where a is the riskless interest rate of return on the firm per unit time, C is 
the total amount payout by the firm per unit time to either the shareholders 
or claims-holders (e.g. dividends or interest payments) if positive, and it is 
the net amount received by the firm from new financing if negative. The 
constant L represents the past length while T is the maturity date. The 
function 5^ : M ^ M is a continuous representing the volatility function on 
the firm value per unit time. The initial process Lp : Vt ^ (7([— L, 0],]R) is 
J^o-measurable with respect to the Borel a-algebra of (7([— L, 0], M), actually 
if is the past value of the firm. The process Vl^ is a one dimensional standard 
Brownian motion adapted to the filtration (^t)o<t<T- 

Notice that C and a can be time dependent functions, in which case they 
should be measurable and integrable in the interval [0,T]. 

The results ensuring the feasiblility of the price model (fTl) is given in 
[Tllin]. Following the work in [7J, in order the RPDE which must be satisfied 
by any security whose value can be written as a function of the value of the 
firm and time, we assume that any claim with market value Y{t) (which can 
be replicated using self-financed strategy) at time t with Y{t) = F{V{t)^t) 
follows a nonlinear stochastic delay differential equation 

dY{t) = {ayY{t) - Cy)dt + gy{Y{t - L))Y{t)dWy{t), t G [0, r] 

(2) 
Y{t) = ^y{t), tG[-L,0], 

on a probability space (f], J^, P)). where ay is the constant riskless interest 
rate of return per unit time on this claim; Cy is the amount payout per 
unit time to this claim; 5^^ : M ^ M is a continuous function representing 
the volatility function of the return on this claim per unit time; the initial 
process (py : Q ^ (7([— L,0],]R) is J^o-i^easurable with respect to the Borel 
cr-algebra of C([— L, 0],]R). The functions Cy{t) and ay{t) are measurable 
and integrable in the interval [0,T]. The process Wy is a one dimensional 
standard Brownian motion adapted to the filtration {J~'t)o<t<T- For any claim 
Y{t) = F{V{t)^t) where F is twice continuously differentiable with respect 
to V and once differentiable with respect to t, we have proved in [14] that 



the following (RPDE) should be satisfied 

^g\V{t - L)yF,, + {rv - C)F, + Ft-rF + Cy = 0, (3) 

{t, v)e[T- L, T] X M+ (4) 

with 

By setting 

V{t) = F{V{t),t) + f{V{t),t), (5) 

where f(y{t)^t) is the value of the equity, F(y{t)^t) the value of debt a any 
time t before the maturity and r is the instantaneous riskless rate of interest. 
We have obtained in [14j the following two final value problems for debt and 
equity, linked by ^ 

\g\V{t - L)yF,, + {rv - C)F, + Ft - rF + Cy = 0, t e {T - L,T) 

F{v,T) = mm[V,B], v>0 (6) 

, F(0,t) = 0, F{v,t) - 5e-^(^-^), as ^ ^ oc, 
and 

^g^Vit - L)yU. + {rv - C)U - rf + C - Cy + ft = 0, t e {T - L, T) 

f{v, T) = max{v - B,0), v>0 C^) 

^ /(O, t) = 0, f{v, t)^v- 5e-'^(^-*), as w ^ oo, 

where B is the promised value the firm must pay to the debtholders at the 
maturity date T. 

Remark 2.1. The classical Merton model I? J assumes that the value of the 
firm at time t G [0, T], V{t), follows a Stochastic Differential Equation (SDE) 

dV = {aV - C)dt + aVdW 

(8) 
y(0) given, 

6 



where a^ is the constant instantaneous variance of the return on the firm 
per unit time. In this case, the equity value f should satisfy the following 
deterministc PDF 



Icr'v'fvv + {rv -C)f,-rf + C-Cy + ft = 0,te (0,r) 
f{v, T) = max(^ -B,0),v>0 (9) 

, /(O, t) = 0, f{v, t)^v- Be-^^^-^\ as v^oo. 

3. Numerical techniques and applications 

3.1. Presentation of the data set and volatility estimation 

The data on stock returns come from the Center for Research in Securty 
Prices(CRSP) database: htt p://www.crsp. com/^ while those on debt values 
are from the Research Insight/Compustat database (lhttp://www. compustat.com/) 



More data include firms that had valid data for all 20 years from 1991-2010 
and including: 

1. The risk free rate r per year, which is the average monthly yield on US 
T-Bills for that year (the same for all firms each year). 

2. The standard deviation of daily returns a per year for each firm. 

3. The number of daily returns A^ used to compute a for each firm each 
year (this is set to be at least 150). 

4. The total book value of debt B (in 1,000,000s). 

5. The total value of the firms assets V (in 1,000,000s). 

6. The total amount C (in 1,000,000s) payout by the firm per unit time 
to either the shareholders or claims-holders for 10 years (2000-2010). 

7. The total amount Cy (in 1,000,000s) payout per unit time for the debt 
within 10 years (2000-2010). 



In fact the data set we have used include aU the parameters that we need to 
solve either the stochastic differential equations ([l]) & ([s]), or the RPDE ^ 
k andthePDE Q. 

All the simulation is performed in Matlab 7.7. In most of our simulations, 
the data between 1991-2000.5 are used as memory data while those between 
2000.5-2010 are used as the future data i.e. the data that we want our model 
to predict. 

To estimate the volatility function g^ we use the quadratic or linear inter- 
polation of the memory part of data a. As in [3j, the quadratic form of the 
volatility is motivated by the fact that the implied volatility in Black-Scholes 
model has a parabolic shape. The volatility function g can also be estimateed 
by using the splines interpolation of the memory part of the data a. 

As we only have yearly data set, we use also the interpolation to have 
more data set if need as the numerical schemes usually need small time step 
(then more data set) to ensure their stabilities. 

3.2. Numerical approximation of the corporate market value 
3.2.1. The 9— semi implicit Euler-Maruyama scheme 

Here we consider the stochastic equations ([l]) and ([s]) within the time 
interval [0,T], where the higher value of T is 9.5 corresponding to the year 
2010. The time unit being the year. Indeed the values of a and C are time 
depending, we therefore consider those values as two time depending func- 
tions, which are constant within each year interval. The goal here is to use 
the mean numerical solutions of the firm V as the forecasting values of firm in 
the interval [2000.5, 2000.5 + T]. As the real firm value of the companies are 
already known in that interval, the aim is to see how close are the forecast- 
ing firm values (from numerical methods) comparing to the real firm values 
from financial industries. Recall that our nonlinear model used the memory 
data within the interval [1991,2000.5]. We solve numerically our nonlinear 
model ([l]) for the value of the company and Merton model ([s]) [j in time 
using the 6— implicit Euler Maruyama scheme in order to obtain numerical 
sample solutions. Monte carlo method is thereafter used to approximate the 
mean numerical solution of each model. The mean numerical value from each 
model will be therefore compared with the real company value V . 



^ For constant a and C the exact solution is well known as this is the same as Black 
Sholes model for stock price 



The 9— semi implicit Euler-Maruyama scheme apphed to ([T]) is given by 

K+i = K + Ar[e(a,+iK+iK-m+i-a+i) + (i-^)KKK-m-C'n)] 

+g{Vn-m)VnAWn n = l , M, 0<9<1, L = mAr. (10) 

where AT = T/M is the time step size, M the total number of time subdivi- 
sion, Vn is the approximation of V{tn)^ tn = nAT, an = (^{tn)^ Cn = C{tn)^ 
and 

are standard Brownian increments, independent identically distributed 
a/ATA/'(0, 1) random variables. For ^ = 0, we have the classical Euler- 
Maruyama scheme which is less numerical stable than the semi implicit 
Euler-Maruyama with = 1, that we will use in our simulation. To en- 
sure the convergence of the numerical (10) toward the unique solution of ([l]), 
the volatility function g need to be globally Lipschitz, or localy Lipschitz 
and bounded [16j. These conditions are sufficient conditions for the conver- 
gence and not necessary conditions since the scheme can converge for some 
functions not verifying these conditions. 

To approximate the expected value (mean) of the process V ^ we use the 
Monte Carlo method to compute the mean of the numerical samples from 



(10). The Monte Carlo method can also be used to approximate any moment 



of the process V . 

3.2.2. Application with corporate data 
The following firms are used: 



Ci Great Northern Iron Ore Pptys (Figure 3(a) and Figure |3(b)D 

C2 Tor Minerals Intl Inc (Figure [3(c)] and Figure |3(d)D 

C^ South Jersey Inds Inc (Figure 2(a) and Figure [2(b)" ) 

(^4 Rentech Inc (Figure 1(c) and Figure [T(d)] ) 

Cs Magna International Inc (Figure [l (a) [ and Figure [T(b)] ) 

Cq First Citizens Bancshares Inc NC (Figure 2(c) [ and Figure [2(d)] ) 



As we have already mentioned, the time origin corresponds to the year 
(2000+1/2), the data before are memory data and we want to predict the 
data after (2000+1/2). We plot 400 samples of the numerical solution for 
our delayed model and Merton model along with the means of the numerical 
samples (green curves). As the origin is year (2000+1/2), the part of the 
mean curves before the origin are just the curves of the real firm market 
value V in that interval. The curves of the real firm market value 1/ as a 
function of time are in black (black thick curves) . Indeed we want the means 
of the numerical samples (green curves) to fit well the real firm market value 
V (black thick curves) with moderate standard derivations (few spread of 
the numerical samples comparing to its mean), this will be the aim of our 
comparisons. 

In all graphs, the function g (volatility in delayed model) is the quadratic 
interpolation of the standard deviation of daily returns a in the memory part 
while the volatility in the Merton model is just the mean of the memory part. 

In Figure [T] we take L = T = 9.5, the graphs at the left hand size. 
(Figure 1(a) and Figure 1(c) respectively for firms C^ and C4 ) correspond to 



the delayed model while the graphs at the right hand size (Figure 1(b) and 



Figure |l(d) respectively for firms C^ and C4) correspond to Merton model. 
For corporate Cs, Figure |l(a)| shows the good prediction with reasonable 
standard deviation (as the numerical samples are not much spread) of the 



delay model while Figure [T(b) shows the early good prediction of the Merton 
model but the prediction has f allied just after the year 2005. For corporate 
C4, Figure 1(b) [ and Figure 1(d) [ show the good prediction before 2005 with 
relatively large standard deviation of the delayed model and Merton model. 

In Figure [2} we aslo take L = T = 9.5, the graphs at the left hand size 
(Figure [2 (a) I and Figure 2(c) [ respectively for firms Cs and Cq) correspond to 
the delayed model while the graphs at the right hand size (Figure |2(b) and 
Figure [2(d)] respectively for firms C3 and Cq)) correspond to Merton model. 
For corporate Ca, Figure [2(a)] shows an early accepted prediction for delayed 
model comparing to the Merton model in Figure [2(b)[ where the prediction 
is more bad (black thick curve been really far away from green curve). For 
corporate Ce, we have the same observation as for corporate Cs according to 
Figure 2(c) and Figure [2(5)" . 

In Figure [3] we take L = 9.5, T = 5, the graphs at the left (Figure 3(a) 



and Figure 3(c) respectively for firms Ci and C2) correspond to the delayed 



model while the graphs at the right (Figure 3(b) [ and Figure [3(d)] respectively 
for firms Ci and C2) correspond to Merton model. From Figure [3l we can 
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observe that for corporate Ci the delayed model fit well the real data of the 
firm market value compared to the Merton model, while for corporate C2 the 
two models fit well the real data of the firm market value before 2003. 

3.3. Numerical Evaluation of Debt or Equity in a Levered Firm 
3.3.1. Numerical schemes based on exponential integrators 

Debt and equity are linked by relation ([5]) , so we only need to solve one of 
the systems ([T]) and ^ . We consider here the random partial differential ([T]) , 
but where C, Cy and r are time depending functions. In our simulation we 
consider those values as time depending functions, which are constant within 
each year interval as we have in our data set. Indeed to solve numerically 
this equation the domain of v need to be troncated. Taking in to account 
the fact that (7, Cy and r are time depending functions, we therefore have 

i g\V{t - L)y /,, + {r{t)v - C{t))U + ft- r{t)f - Cy{t) + C{t) = 0, 

f{v,T) = max{v - B,0), v G [0, K,ax] (H) 

/(0,t)=0, te[T-L,T] 

f{Vra...t) = V^,^-Be--^tr{s)ds^ t ^ [T - L,T] 



Our model problem (11) is similar to Europeans call options prices, we can 
therefore take Knax three or four times B according to \Xh\' In our simulation 
we take 14iax = 45, where B is the amount that the firm must pay to the 
debtholders at the maturity date T (like the strike price for options prices). 



Our system (11) is a backward system, to transform it to the forward one, 
we use the tranformation r = T — t^ and the corresponding equation is given 

by 

i g\V{T -T- L)y U + {r{T)v - C{t))U - r{r)f 

-Cy{T)+C{T)=fr, 

f{v, 0) = max(t- - 5, 0), vE [0, KnaJ (^2) 

/(0,r) = 0, re[0,L] 

Please note that after the transformation r = T — t^ the functions C(r), C^(r) 
and r{T) in (12) are normally the functions C{T — T)^Cy{T — r) and r{T — r). 



To apply sophistical technique to the convection term (the term with /^) 
in order to avoid numerical instabilities, let us put this term in the so called 

11 



the conservation form. In fact 

{r{T)v - C{t))U = ir{r)v - C(r))/), - r{T)f 
Using this relation, equation ( [T2| ) become 

^ I g\V{T - T - L)y U + {{r{r)v - C{r))fX - 2r{r)f, 

-CyiT) + CiT) = f^ 

/(0,r)=0, re[0,L] 



To solve equation ([13|) two cases can be considered: 

1. The case where T — r — L < 0, Vr G [0, L], then T < L. 

2. The case where T > L. 



For the first case (T < L) the RPDE (13) become the deterministic PDE 
since V{t) = Lp{t) for t G [— L, 0] as given in ([l]). 

For the second case {T > L)^ to solve (13) the following step should be 
followed 

1. Solve the stochastic equation ^ to have a sample of the numerical 
solution of V as we did in the previous section. 



2. Use the numerical sample solution of V from step 1 to build the diffu- 



sion coefficient (the coefficient of fyy) in the RPDE (13), which therefore 



become a deterministic PDE for this fixed numerical sample of V. 

3. Solve the deterministic PDE from step 2 for the fixed numerical sample 
of V from step 1. 

4. Repeat step 1, step 2 and step 3, M times (relatively large) and use 
the Monte Carlo technique to estimate the expectation value of / and 
also any moment of the stochastic process / if need. 
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As the two cases require the solution of the deterministic PDE, in the sequel 
we will consider the first case (T < L)^ and the corresponding deterministic 
PDE is given by 

^ I g\^{T - T - L)y /,„ + {{r{r)v - C{r))f)^ - 2r{r)f 

-Cy{T)+C{T)=U, 

f{v,0) = ms.x{v-B,0), ve[0,V^a.] (14) 

/(0,r) = 0, re[0,L] 

For the discretization in the direction of 'U, we use the combined finite 
difference finite volume method. The interval [0, Vmax] is subdivised into 
A^ parts that we assume equal without loss the generality. As in center finite 
volume method, we approximate / at the center of each interval. The diffu- 
sion part of the equation is approximated using the finite difference while the 
convection term is approximated using the standard upwinding usual used 
in porous media fiow problems [T9l [II [Ml |20] . 
Let 

v, = {2i-l)h/2, h=^, z = l,2, ,A 

being the center of each subdivision . We approximate the diffusion term at 
each center by 

lg\^{T-r-L))v^fUv^) ^ ^^ g\^{T - r - LM x 

ifi+^ir) - 2Mt) + fi.,{T)) , ^ = 2, ..., iV - 1. 
\g\^{T - T - L))vl UM) ^ ^g'{^{T-r-L))vlx 

/2(r)-/i(r) Mr) 
h h 

^g^ip{T-T-L))v%U{vM) ^ ^g\ip{T-T-L))v%x 

f f{Vmax,r) - fNJr) _ fNJT) - /jv-i(r) 

V h/2 h 

This approximation is similar to the one in [21j with central difference on 
non uniform grid. We approximate the convection term using the standard 
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upwinding technique as following 

. . , rf x,n ( ^^ (K^)^.+i/2 - C(r))/+(r) - {riT)v,-i/2 - C{T))fUT) 
{r{T)v - C{T))f)^ [Vi) ^ 

where 

f+(^\ _ S fii^) if r{T)vi+i/2-C{T)^0 , . 

^'^^ ~ Xfi-iir) if r(rK+i/2 - C(r) < ^'"^ 

Vi+i/2 = Vi + h/2, Vi_i/2 = Vi- h/2 = Vi-i + h/2, (16) 

where 

Mr) ^ fivi,T). 

Reorganizing all previous diffusion and convection approximations lead 
to the following initial value problem 

^ = A(r)f + b(r), re[0,L] ^^^^ 

f (0) = (max('Ui — B, 0), ..., max('UAr — B, 0)) . 

where A(r) is a tridiagonal matrix and 

f (r) = (/.(r))i<,<^ , b(r) = C{t) - Cyir) + k(r). (18) 

where k is the contribution from boundary conditions. 

The function x h^ max(x, 0) is not smooth, it important to approximate 
it by a smooth function. The approximation in [2T] is a fourth-order smooth 
function denoted tt^ and defined by 

{X if X ^ e 
Co + Cix + + Cgx^ if — e < X < e (19) 

if x<-e 

where < e <C 1 is the transition parameter and 

35 1 35 35 

<^0 = 7^^^, Ci = -, C2 = -7—, C4 



256 ' ' 2' ^ 64e' " 12863' 

This approximation allow us to write 

f(0) = 7r,(v-i?), v=(i;,),<,<^. (20) 
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Let us introduce the time stepping discretization for the ODE (17) based 



on exponential integrators. Classical numerical methods usually used are 
Implicit Euler scheme and Crank-Nicolson scheme ^22j. Following works 



from [211 [ini U] the exact solution of (17) is given by 

rrn+Ar 






Tn + Ar 



^JTn 



A(s)ds 



J Tr, 



(21) 



r^ = nAr, n = 0, ...,M, Ar > 0. (22) 



Note that (21) is the exact representation of the solution, to have the nu- 



merical schemes, approximations are needed, the first approximations (using 
the quadrature rule) may be 

/ A(s)ds ^ AtA(t„) / A{y)dy^{s-Tn)A{Tn) (23) 

Using these approximations we therefore have the following second-order 
approximation 



f (r. + Ar) 



^ArA(rn) 



rrn+Ar 



f (r^) + / ^ e-^'-^-^^^^-^h{s)ds 



(24) 



The simple scheme caUed Exponential Differential scheme of order 1 (ETDl) 
is obtained by approximating b(s) by the constant b(r„) and is given by 



in+i = fn + (ArA(r„))-' [A(r„)f„ + b(r„)] . 



(25) 



A second order scheme is given in [21j. 

Following the work in [23^^ Lemma 4.1], if the the function b can be well 
approximated by the polynomial of degree p (which is the case here since we 



have the exponential decay at the boundary v = Vmax)^ from (24) we have 



in+i = M^TA{Tn)% + J]^-^Ar'+V;+i(ATA(T„))b,+i, (26) 



j=0 1=0 



u-iy- 



where 



p-i 



^^^^ ^ I^TT'^^+I' 



ifo{x) 



j=0' 



e"", ifi{x)=xifi+i{x) + -, 1 = 0,1,2, 
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Note that to have high order accuracy in time for p > 2, the integral in (23) 
should be approximated more accurately. The Magnus expansion may also 
used in such case. 

3.3.2. Application with corporate data 

All schemes here can be implemented using Krylov subspace technique in 
the computation the expomential functions presented in those schemes with 
the Matlab functions expmvp.m or phipm.m from j23l [24] . The Krylov sub- 
space dimension we use is m = 10 and the tolerance using in the computation 
of the expomential functions 99^ is to/ = le — 6. We use p = 2 and obtain 



second order accuracy in time as the approximations (23) are second order 
in time. 

We used the following frims: 



C^ Magna International Inc (Figure]^ 

Cq First Citizens Bancshares Inc NC (Figure [^ 

C7 Coca-Cola CO (Figure |6|) 

Cs One Liberty Properties INC (Figure [?]) 

Cq Cisco Systems INC (Figure [^ 

Cio C B S Corp NEW (Figure |9 



Cii Nam Tai Electronics INC (Figure 10) 



Here again, the time origin corresponds to the year (2000+1/2), the data 
before are memory data and we want to predict the data after (2000+1/2). 
In the legends of all of our graphs we use the following notation 

• "Delayed Equity" is for the numerical equity value from our nonlinear 
delayed model. 

• "Real Equity" is for the real equity value of the corporate. 

• "Merton Equity" is for the numerical equity value from Merton model. 
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In our surface graphs of the numerical equity value, we plot only the part 
where the variable V is between the minimun and the maximum values of 
our real market value V. In all simulations with our delayed model, we take 
L — 9.5. In all graphs, the function g (volatility in delayed model) is the 
quadratic interpolation of the standard deviation of daily returns a in the 
memory part while the volatility in the Merton model is just the mean of the 
memory part. 

For each firm, we plot at the left hand size both the surface graphs of the 
numerical equity value from our delayed model at T = 9.5 and T = 5. In 
those 3D surface graphs, we also plot the corresponding 3 D graphs (green 
curves) of the real data of the firm equity value as a function of the time 
(year) and V. At the right hand size, we plot in 2 D the firm equity value 
as a function of time (year) , corresponding to the surface graphs at the left 
hand size. Those 2 D equity graphs contain the numerical equity value from 
our delayed model, the numerical equity value of the Merton model and the 
real data equity value of the firm. 

In our simulations, for a given T, the promised debt B is just the real 
debt value of the firm at time T. 

For firm C^ in Figure [4| we can observe that both the delayed model and 
Merton model fit well the real market equity value of the firm. The accuracy 
of the two methods varies within some time interval as we can observe in 
Figure [4(b)] and Figure [4(d)| 

For firm Cq in Figure [5} comparing to firm C^ the two models fit less. 
In a wide time interval in Figure [5(b)] and Figure 5(d), the delayed model is 



more close to the real market equity of the firm. We can also observe a good 
early fit in the Merton model. 

For firm C7 in Figure [6| comparing to firm C^ the two models fit less. 
But for the maturity date T = 9.5 in Figure [6] the fitting is relatively good 
for the two models. The accuracy of the two methods varies within some 
time interval as we can observe in Figure [6(b)[ and Figure [6(d)[ 

For firm Cs in Figure [7| the fitting is relatively bad for the two models. 
However we can observe in Figure [7(b)] and Figure [7(d)] the good early fit in 
the Merton model, and that in the wide time interval the delayed model is 
more close to the real market equity of the firm than the Merton model. 

For firm Cg in Figure [8| the fitting is relatively good for the two models 
in the early time interval and become relatively bad just after. 

For firm Cio in Figure [9j the fitting is relatively good for the two models 
for the maturity date T = 9.5 in Figure [9(b) at the middle time interval and 
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bad for the maturity date T = 5 in Figure 9(d)[ 



For firm Cu in Figure [TOl the fitting is relatively good for the two models 
in the early time interval but become bad just after. The two models are 
confused. 

4. Conclusion 

In this paper, numerical techniques to solve delayed nonlinear model for 
pricing corporate liabilities are provided. The numerical technique to solve 
the RPDEs modeling debt and equity value combines the finite diflFerence 
finite volume methods (discretization respect to the firm value V) and an 
exponential integrator (discretization respect to the time t). The matrix 
exponential functions are computed efficiently using Krylov subspace tech- 
nique. 

Using financial data from some firms, we compare numerical solutions 
from both our nonlinear model and classical Merton to the real firm's data. 
This comparaison shows that our nonlinear model behaves very well. We 
conclude that in corporate finance the past dependence of the firm value 
process may be an important feature and therefore should not be ignored. 
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2000 
Year 



(a) 



(b) 





(c) 



(d) 



Figure 1: The graphs at the left ((a),(c) respectively for corporates C3 and Cq ) correspond 
to the delayed model while the graphs at the right ((b), (d) respectively for corporates C3 
and Cq) correspond to Merton model. We aslso take T = L = 9.5 and the function g is 
the quadratic interpolation of the standard deviation of daily returns a in the memory 
part. We have plotted 400 samples of the numerical solution along with the expectation 
(the means) of the numerical solution (green curves). The curves of the real data of the 
firm market value F as a function of time are in black (black thick curves). 
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(a) 



(b) 





(c) 



(d) 



Figure 2: The graphs at the left ((a),(c) respectively for corporates C5 and C4 ) correspond 
to the delayed model while the graphs at the right ((b), (d) respectively for corporates C5 
and C4) correspond to Merton model. We also take T = L = 9.5 and the function g is the 
quadratic interpolation of the standard deviation of daily returns a in the memory part. 
We have plotted 400 samples of the numerical solution along with the expectation (the 
means) of the numerical solution (green curves) . The curves of the real data of the firm 
market value F as a function of time are in black (black thick curves) 
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Year 
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Year 



(a) 



(b) 
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Year 




1992 1994 1996 1998 2000 2002 2004 



(c) 



(d) 



Figure 3: The graphs at the left ((a) and (c) respectively for corporates Ci and C2) 
correspond to the delay model while the graphs at the right ((b) and (d) respectively for 
corporates Ci and C2) correspond to Merton model. We take T = 5 L = 9.5 and the 
function g is the quadratic interpolation of the standard deviation of daily returns a in 
the memory part. We have plotted 400 samples of the numerical solution along with the 
expectation (the means) of the numerical solution (green curves). The curves of the real 
data of the firm value V as a function of time are in black (black thick curves). 
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Figure 4: The graphs for Firm C5. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 5: The graphs for Firm C^. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 6: The graphs for Firm C7. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 7: The graphs for Firm C^. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 8: The graphs for Firm Cg. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 9: The graphs for Firm Cio- We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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Figure 10: The graphs for Firm Cu. We plot at the left ((a) and (c) respectively) the 
surface graphs of the numerical equity from our delayed model at T = 9.5 and T = 5. 
The corresponding 3 D graphs (green curves) of the real data of the firm equity value as 
a function of the time and V are also plotted in (a) and (c). At the right ((b) and (d)), 
2 D graphs of the firm equity value as a function of time, corresponding to the surface 
graphs at the left ((a) and (c)) respectively) are presented. Those 2 D graphs contain the 
numerical equity from our delayed model, the numerical equity of the Merton model and 
the real data equity of the firm. 
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